Sparse Convolution Quadrature for Time Domain Boundary Integral Formulations of the Wave Equation
نویسندگان
چکیده
Many important physical applications are governed by the wave equation. The formulation as time domain boundary integral equations involves retarded potentials. For the numerical solution of this problem we employ the convolution quadrature method for the discretization in time and the Galerkin boundary element method for the space discretization. We introduce a simple a-priori cutoff strategy where small entries of the system matrix are replaced by zero. The threshold for the cutoff is determined by an a-priori analysis which will be developed in this paper. This method reduces the computational complexity for solving time domain integral equations from O ( M2N logN ) to O ( M1+sN logN ) for some s ∈ [0, 1[, where N denotes the number of time steps and M is the dimension of the boundary element space.
منابع مشابه
Rapid Solution of the Wave Equation in Unbounded Domains: Abridged Version
Abstract We propose and analyze a new fast method for the numerical solution of time-domain boundary integral formulations of the wave equation. Discretization in time is achieved by Lubich’s convolution quadrature method and in space by a Galerkin boundary element method. We show that the arising block Toeplitz system is after a small perturbation equivalent to a a decoupled system of discreti...
متن کاملNumerical treatment of Retarded Boundary Integral Equations by Sparse Panel Clustering
We consider the wave equation in a boundary integral formulation. The discretization in time is done by using convolution quadrature techniques and a Galerkin boundary element method for the spatial discretization. In a previous paper, we have introduced a sparse approximation of the system matrix by cutoff, in order to reduce the storage costs. In this paper, we extend this approach by introdu...
متن کاملNumerical resolution of large deflections in cantilever beams by Bernstein spectral method and a convolution quadrature.
The mathematical modeling of the large deflections for the cantilever beams leads to a nonlinear differential equation with the mixed boundary conditions. Different numerical methods have been implemented by various authors for such problems. In this paper, two novel numerical techniques are investigated for the numerical simulation of the problem. The first is based on a spectral method utiliz...
متن کاملFast convolution quadrature for the wave equation in three dimensions
This work addresses the numerical solution of time-domain boundary integral equations arising from acoustic and electromagnetic scattering in three dimensions. The semidiscretization of the time-domain boundary integral equations by Runge-Kutta convolution quadrature leads to a lower triangular Toeplitz system of size N . This system can be solved recursively in an almost linear time (O(N logN)...
متن کاملSpace - time BIE methods for non homogeneous exterior wave equation problems . The Dirichlet case . ∗
In this paper we consider the (2D and 3D) exterior problem for the non homogeneous wave equation, with a Dirichlet boundary condition and non homogeneous initial conditions. First we derive two alternative boundary integral equation formulations to solve the problem. Then we propose a numerical approach for the computation of the extra “volume” integrals generated by the initial data. Finally, ...
متن کامل